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We report on the results of multi-epoch VLBI observations with VERA (VLBI Exploration of Radio 
Astrometry) of the 22 GHz H2O masers associated with the young stellar object SVS 13 in the NGC 1333 
region. We have carried out phase-referencing VLBI astrometry and measured an annual parallax of the 
maser features in SVS 13 of 4.25±0.32 mas, corresponding to the distance of 235±18 pc from the Sun. 
Our result is consistent with a photometric distance of 220 pc previously reported. Even though the maser 
features were detectable only for 6 months, the present results provide the distance to NGC 1333 with 
much higher accuracy than photometric methods. The absolute positions and proper motions have been 
derived, revealing that the H2O masers with the LSR (local standard of rest) velocities of 7-8 km s" 1 are 
most likely associated with VLA4A, which is a radio counterpart of SVS 13. The origin of the observed 
proper motions of the maser features are currently difficult to attribute to either the jet or the rotating 
circumstellar disk associated with VLA4A, which should be investigated through future high-resolution 
astrometric observations of VLA4A and other radio sources in NGC 1333. 

Key words: Astrometry: — ISM: individual (NGC 1333) — ISM: jets and outflows — masers (H2O) 
— stars: individual (SVS 13) 



1. Introduction 

The reflection nebula NGC 1333 is one of the near- 
est low-mass star-forming regions and is associated with 
the classical Hcrbig-Haro (HH) objects HH 7-11 (Hcrbig 
1974; Strom ct al. 1974). Despite the importance 
of NGC 1333 for studying star-formation processes and 



the overall structure of the larger molecular cloud com- 
plex, the distance to NGC 1333 is still quite uncertain. 
NGC 1333 is located in the Perseus complex of dark 
clouds, which forms a chain of several clouds with a 
size of 7° x 3°, elongated perpendicular to the Galactic 
plane (Cernis 1990, 1993). The molecular clouds asso- 
ciated with NGC 1333 are located at the western edge 
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of the chain, while the young open clusters IC 348 and 
Perseus OB2 are at the eastern edge of this chain (Cernis 
1990, 1993). 

The large scale structure of the Perseus dark cloud com- 
plex has been studied extensively (Cernis 1990, 1993; de 
Zeeuw et al. 1999). According to these studies, there 
is a gradient in distance across the complex (see Fig. 9 of 
Cernis 1993); NGC 1333 is the nearest cloud at a distance 
of 220 pc (Cernis 1990) and IC 348 and the Perseus OB2 
association are at a distance of 318-340 pc (Cernis 1993; 
dc Zceuw et al. 1999). Although Hipparcos provided an 
accurate distance to the Perseus OB2 association, based 
on an annual parallax measurement of 318±27 pc (de 
Zeeuw et al. 1999), the distance to NGC 1333 and the 
other dark clouds have been measured only through opti- 
cal photometry, with a typical uncertainties of about 25% 
(Cernis 1990, 1993). Therefore, further precise astromc- 
try allowing an accurate annual parallax measurement is 
necessary to reveal overall structure of the chain of the 
molecular clouds in the Perseus region. 

The young stellar object (YSO) SVS 13 is located at the 
base of HH 7-11 in NGC 1333 and has been proposed to 
be a powering source of the jets and outflows in this region 
(Strom et al. 1976). High-resolution observations with 
radio interferometers at centimeter and millimeter wave- 
lengths have been conducted to investigate the complex 
nature of this system. Very Large Array (VLA) obser- 
vations at 3.6 cm wavelength reveal a radio continuum 
source VLA4 that appears to be associated with SVS 13 
(Rodriguez et al. 1997, 1999). Higher resolution VLA 
observations of SVS 13 at the same wavelength resolved 
VLA4 into a double radio source, denoted VLA4A and 
VLA4B (Anglada et al. 2000), with the separation angle 
of 0.3". Anglada et al. (2000) reported that the position 
of SVS 13 coincides with the western component, VLA4A, 
while that of the millimeter source (Looney et al. 2000) is 
associated with the eastern component, VLA4B. VLA4B 
is suggested to have larger amount of circumstellar mate- 
rial than VLA4A, based on the subsequent high-resolution 
VLA observations in the 7 mm band (Anglada et al. 
2004). 

Although HH 7-11 have been proposed to be powered 
by SVS 13, Rodriguez et al. (1997) argued that another 
radio continuum source VLA3, located 6" southwest of 
SVS 13, is a more favorable candidate for the powering 
source due to the better alignment of VLA3 with the 
HH objects. In contrast, Bachiller et al. (2000) pre- 
sented the results of the intcrferometric CO J=2-l obser- 
vations of HH 7-11 that indicate that the extremely high 
velocity (EHV) molecular outflow is powered by SVS 13, 
rather than VLA3 which has no high velocity CO outflow. 
Looney et al. (2000) also suggested that VLA3 (or mil- 
limeter source A2 in their paper) is a less likely candidate 
for the source that powers HH 7-11. Rodriguez et al. 
(2002) claimed that both VLA3 and SVS 13 (VLA4B is 
favorable rather than 4A) could be the powering source 
of HH 7-11. Despite the fact that the axes of the HH jets 
and molecular outflows are similar, their exciting sources 
could be different, since there are several YSOs in a small 



area. In fact, it is possible that there is more than one 
molecular outflow because, as noted by Rodriguez et al. 
(2002), inspection of Figure 2 of Bachiller et al. (2000) 
shows that the EHV outflow is associated with VLA4, 
while the standard high velocity (SHV) gas is found in 
the vicinity of both VLA3 and VLA4. Thus the origin 
of the jets, molecular outflows, and HH 7-11 is still am- 
biguous and might be revealed through accurate proper 
motion measurements. 

Proper motions can be accurately measured by phasc- 
rcfcrcncing VLBI observations. If one employs cxtragalac- 
tic radio sources as the position references, one can mea- 
sure the absolute position of the target source (Beasley & 
Conway 1995), making it feasible to derive its annual par- 
allax as well as the absolute proper motion. Highly pre- 
cise VLBI astrometric observations have been carried out 
for the Galactic CH3OH and H2O maser source W3(OH) 
with the NRAO Very Long Baseline Array (VLB A), which 
yielded distances of 2.0 kpc from the Sun with uncertain- 
ties of about 2%, based on the annual parallax method 
(Xu et al. 2006; Hachisuka et al. 2006). 

Recently, we have constructed a new VLBI network 
in Japan called VERA, VLBI Exploration of Radio 
Astrometry (Kobayashi et al. 2003), which is the first 
VLBI array dedicated to phase-referencing astrometry. 
Initial results of VERA observations have been reported 
(e.g. Sato et al. 2007; Honma et al. 2007; Hirota et 
al. 2007; Imai et al. 2007), in which absolute proper 
motions and annual parallaxes are measured for galactic 
H2O maser sources at distances ranging from 178 pc for 
IRAS 16293-2422 (Imai et al. 2007) to 5.28 kpc for S269 
(Honma et al. 2007), demonstrating VERA's astrometric 
capability. 

In this paper, we present the results of astrometry of 
the H2O masers associated with SVS 13 in NGC 1333 
with VERA. These observations have been done as one of 
the initial scientific projects of VERA: measurements of 
annual parallaxes of nearby molecular clouds (e.g. Hirota 
et al. 2007; Imai et al. 2007). SVS 13 is known to be one 
of the brightest H2O maser sources among the known low- 
mass YSOs (Claussen et al. 1996; Furuya et al. 2003). 
In addition, there is a bright ICRF source, J0336+3218, 
which is separated from SVS 13 by an angle of 1.89° (Fey 
et al. 2004). The aim of this paper is to measure the 
distance to NGC 1333 as well as to explore the powering 
source of the jets and outflows in this region through the 
maser positions and proper motions. 

2. Observations 

The VERA observations of H2O masers (6i6-523, 
22235.080 MHz) associated with SVS 13 were carried out 
for 10 hours in each of 7 observing sessions from Nov. 2004 
to May 2005 (2004/317, 2004/355, 2005/021, 2005/052, 
2005/080, 2005/113, and 2005/140; hereafter an observ- 
ing session is denoted by year/day of the year) at intervals 
of about 1 month. All 4 stations of VERA (see Fig.l of 
Pctrov et al. 2007) took part in all sessions, providing a 
maximum baseline length of 2270 km. 
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Observations were made in the dual beam mode; the 
H2O masers associated with SVS 13 and a reference source 
J0336+3218 (Fey et al. 2004), with the separation angle 
of 1.89°, were observed simultaneously. The instrumental 
phase difference between the two beams was measured 
continuously during the observations by injecting artificial 
noise sources into both beams at each station (Kawaguchi 
et al. 2000; Honma et al. 2003). The typical value of 
the phase drift between the two beams was 3 degrees per 
hour, which was removed from the data. 

Left-handed circular polarization was received and sam- 
pled with 2-bit quantization and filtered using the VERA 
digital filter unit (Iguchi et al. 2005). The data were 
recorded onto magnetic tapes at a rate of 128 Mbps, 
with two IF channels of 16 MHz bandwidth each for both 
SVS 13 and J0336+3218. J0336+3218 was detected with 
a peak intensity greater than 1.2 Jy beam -1 in all the ses- 
sions, so it was also used for bandpass and delay calibra- 
tion. System temperatures and atmospheric attenuation 
were measured with the chopper-wheel method (Ulich & 
Haas 1976), and typical values were 100-400 K, depending 
on weather conditions and elevation angles of the observed 
sources. The aperture efficiencies of the antennas ranged 
from 45 to 52%. Correlation processing was carried out on 
the Mitaka FX correlator (Chikada ct al. 1991) located 
at the NAOJ Mitaka campus. For H2O maser lines, the 
spectral resolution was set to be 15.625 kHz, correspond- 
ing to the velocity resolution of 0.21 km s -1 . 

3. Data Reduction 

Data reduction was performed using the NRAO 
Astronomical Image Processing System (AIPS). 
Amplitude and bandpass calibration were done for 
the target (SVS 13) and reference (J0336+3218) sources 
independently. For phase calibration, fringe fitting was 
done with the AIPS task FRING on the phase reference 
source (J0336+3218), and the solutions were applied to 
the target source (SVS 13). In addition, we applied the 
results of the dual-beam phase calibration measurements 
(Kawaguchi ct al. 2000). We also corrected for the 
approximate delay model adopted in the correlation 
processing (Honma et al. 2007) and for drifts of the 
visibility phase caused by the Earth's atmosphere based 
on the GPS measurements. 

We first searched for H2O maser "spots", defined as 
emission occurring in a single velocity channel, in fringe- 
rate maps with the AIPS task FRMAP. We found that the 
H 2 masers were detected toward two different sources 
in NGC 1333; one was associated with SVS 13, which 
was previously identified to be H2 0(A) by Haschick et 
al. (1980), and another was located 34" southwest from 
SVS 13, corresponding to the radio continuum source 
VLA2 and H 2 0(B) by Haschick et al. (1980). The H 2 
masers associated with VLA2 were detected only in the 
first three sessions from 2004/317 to 2005/021. Because 
the positions of the maser spots associated with VLA2 
are shifted by 34" southwest of the phase-tracking center, 
the synthesized images are significantly distorted due to 
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Fig. 1. Spectra of the H2O maser features associated with 
SVS 13. Solid and dashed lines represent the features 1 
and 2, respectively. Error bars represent the fitting errors 
(la) given by the AIPS task JMFIT. Maser spots with sig- 
nal to noise ratios larger than 5 times the noise (5<r) for at 
least two consecutive channels are considered to be signif- 
icant detections and are plotted with open circles for fea- 
ture 1 and filled squares for feature 2 (see text). For lower 
signal to noise (non-detections) we only plot error bars. 
Note that negative intensities are expected due to thermal 
noise variations. Spectra at different epochs are offset by 
5 Jy beam -1 steps to clarify the time variations of the spectra. 
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the time-averaging smearing effect (Cotton 1999), mak- 
ing it almost impossible to conduct accurate astrometry. 
Therefore, we do not discuss the masers associated with 
VLA2 in this paper. 

High resolution synthesis imaging and deconvolution 
(CLEAN) were done for the maser spots associated with 
SVS 13 using the AIPS task IMAGR. The image size of 
each map was 51.2 mas x 5 1.2 mas with the pixel size of 
0.1 mas x 0.1 mas. The naturally- weighted synthesized- 
beam size (FWMH) was typically 1.3 masx0.9 mas with 
a position angle of —50°. With a net integration time 
of 7 hours for the H2O maser lines, the resultant rms 
noise levels of the phase-referenced images ranged from 
0.06 Jy beam -1 to 0.16 Jy beam -1 in a single spectral 
channel. We regard the H2O maser spot to be real if the 
signal to noise ratio of the peak intensity is larger than 

5 times the noise level (5a) for at least two consecutive 
channels, provided the positions are coincident within the 
synthesized beam size. In order to improve the signal to 
noise ratios in the phase-referenced images, we made inte- 
grated intensity maps of the maser features, by summing 
the channel maps over the detected velocity range. The 
peak positions and peak intensities of maser features were 
derived by fitting elliptical Gaussian brightness distribu- 
tions to these integrated intensity maps using the AIPS 
task JMFIT. The formal uncertainties in the feature po- 
sitions given by JMFIT were 0.03-0.1 mas, depending on 
the signal to noise ratios of the features and possibly their 
spatial structure. The uncertainties in the peak intensities 
were 0.1-0.3 Jy beam -1 (la) which were approximately 
equal to or slightly larger than the rms noise levels of the 
phase-referenced images. 

4. Results 

4-1. Structure of the H2O maser features 

Figure 1 shows the spectra of the H2O masers associ- 
ated with SVS 13, which are obtained by fitting Gaussian 
brightness distributions to each channel map. We de- 
tected two spatially distinct H2O maser features at the 
LSR velocities of 7-8 km s _1 , which agree well with the 
ambient cloud velocity (e.g. Rodriguez et al. 2002). 
Hereafter we call them as feature 1 and 2, as explained 
later. We did not find other velocity components detected 
previously by Haschick et al. (1980), Wootten et al. 
(2002), and Rodriguez et al. (2002), due to the variabil- 
ity of the H2O masers (e.g. Claussen et al. 1996; Furuya 
et al. 2003). As shown in Figure 2, the peak velocity of 
each maser feature was found to drift systematically from 
7.0 km s -1 to 7.9 km s _1 during the observing period of 

6 months, at a rate of 1.7 km s -1 yr -1 . It is likely that 
the observed velocity drift indicates true acceleration of 
the maser feature as discussed later, although we cannot 
rule out the possibilities of a change in strength of blended 
hyperfine components (Walker 1984) or the structure in 
the maser feature. 

Figure 3 shows the positions of the maser spots and 
features associated with SVS 13. Initially, we found only 
feature 1 in the sessions 2004/317 - 2005/021, as seen in 
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Fig. 2. Peak LSR velocities of the H2O maser features. Open 
circles and filled squares represent the peak velocities for the 
features 1 and 2, respectively, which are derived from the 
weighted means of the spectra shown in Figure 1. Error 
bars indicate the standard deviations (ltr) of the velocities. 



the spectra (Figure 1). In the session 2005/052, feature 
2 appeared at 4.5 mas north of feature 1. The two fea- 
tures existed at almost the same LSR velocities of about 
7.5 km s" 1 in the sessions 2005/052 and 2005/080. In 
2005/080, feature 2 became brighter than feature 1, which 
disappeared in the subsequent session 2005/113. Feature 2 
remained until the session 2005/140, while all the H 2 
masers associated with SVS 13 finally disappeared after 
the 6-month monitoring period of our observations. 

We can see the internal structures of the maser features, 
which are the collections of several maser spots, in Figure 
3. The dispersions of the maser spots within the features, 
0.04-0.27 mas, arc significantly smaller than the synthe- 
sized beam size, except for those detected in the session 
2005/021 in which we found two different groups of spots 
with the separation of 1 mas. These northern spots corre- 
spond to the weak red shifted shoulder at the LSR velocity 
of 7.5-7.7 km s _1 as seen in Figure 1 and can be recognized 
as part of the elongated structure of the maser feature in 
the integrated intensity map. These results suggest that 
the peak position of the maser feature possibly affect the 
astrometric accuracy (e.g. Hirota et al. 2007; Imai et al. 
2007), as discussed in the next section. 

4-2. Astrometry of the H2O maser features 

As shown in Figure 3, the movement of the masers sig- 
nificantly deviates from a simple linear motion; as we will 
show the data can be well modeled by the effects of an 
annual parallax. Assuming that the movement of the 
maser feature is the sum of a linear motion and the annual 
parallax, we derive the proper motions in right ascension 
H a cos 5 and declination /is, the initial positions in right 
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Table 1. Results of the least-squares analysis for the an- 
nual parallax and absolute proper motion measurements 
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Fig. 3. Distribution of H2O maser spots associated with 
SVS 13. Open circles and squares represent the posi- 
tions of the maser spots belonging to the features 1 and 
2, respectively, while crosses represent features that are 
the collections of spots integrated over contiguous spec- 
tral channels. Note that the formal errors given by the 
Gaussian fitting with the AIPS task JMFIT are as small 
as 0.1 mas or less, which cannot be shown in the figure. 
The positions are with respect to the reference position at 
a(J2000) =03h29m03.72465s, <5(J2000) = +31dl6'03.8015". 



ascension ao and declination 80, and the annual paral- 
lax 7r for the maser feature by a least-squares analysis 
as summarized in Table 1 and Figure 4. In the least- 
squares analysis, we determined the proper motions and 
initial positions of two maser features independently, while 
the annual parallax of SVS 13 is common for both fea- 
tures. According to Honma ct al. (2007) and Hirota 
ct al. (2007), the precision of the derived annual paral- 
lax is significantly improved when we use only the data 
for right ascension, as this data is less affected by atmo- 
spheric modeling errors. In this manner, we find the an- 
nual parallax of the H2O maser features associated with 
SVS 13 in NGC 1333 to be 4.25±0.32 mas, corresponding 
to a distance of 235±18 pc. We note that the parallax 
derived from the right ascension data fits the declination 
data quite well. 

According to our astromctric results, the standard devi- 
ations of the post-fit residuals are 0.10 mas and 0.29 mas 
in right ascension and declination, respectively, as listed 
in Table 1. Were we to assume distances of 300 pc and 
350 pc from the Sun (similar to that of IC 348 and the 
Perseus OB2 association), the standard deviation in right 
ascension degrades to 0.19 mas and 0.26 mas, respec- 
tively. Thus, the distance to NGC 1333 is most likely 
235 pc, which is consistent with the photometric distance 
to NGC 1333 reported by Cernis (1990) of 220 pc with an 
uncertainty of 25%, rather than the larger value of about 
300 pc (e.g. Herbig & Jones 1983; de Zeeuw et al. 1999). 



Parameter 



- 1 ) 



a (mas)" 
5q (mas) a 

/2 a cosS (mas yr -1 ) 
fig (mas yr" 
/i (mas yr ] 
Vt (km s _1 ) 
PA (degrees) 
7r (mas) 
o~ a (mas) 
<jg (mas) 



Feature 1 
-4.08(0.14) 



6.14(0.31) 
17.9(0.9) 
-7.9(1.4) 
19.6(1.0) 
21.8(1.1) 
114 
4.25(0.32) 
0.10 
0.29 



Feature 2 
1.17(0.36) 



.38(0.32) 
10.6(1.7) 
10.0(2.1) 
14.6(1.9) 
16.3(2.1) 
133 



Note — Numbers in parenthesis represent the 
estimated uncertainties. Annual parallax tt is 
derived from the right ascension data only. 
a a and 17,5 are the rms deviations 
of the post-fit residuals in the right ascension 
and declination directions. 

a: The position offsets are measured with respect 
to the reference position. 



This is the first time that the distance to NGC 1333 is de- 
termined based on an annual parallax measurement with 
an uncertainty of only 8%. Even though the observing pe- 
riod was as short as 6 months, our results provide a strong 
constraint on the distance to NGC 1333 with higher preci- 
sion than that of the photometric method (Cernis 1990). 

Unfortunately, the masers associated with SVS 13 com- 
pletely disappeared after the 6-month monitoring obser- 
vations with VERA. We will continue long-term monitor- 
ing observations of the H2O maser sources in NGC 1333, 
including SVS 13 and others such as IRAS2 and IRAS4 
(Rodriguez ct al. 2002), which may enable us to improve 
the accuracy of the annual parallax measurements with 
VERA. 

Combining our new results with those of de Zeeuw et 
al. (1999) based on the Hipparcos measurements of the 
distance to the Perseus OB2 association, 318±27 pc, we 
provide a definite evidence for a distance gradient along 
the chain of molecular clouds in the Perseus region as 
proposed by Cernis (1990) and Cernis (1993). Although 
we determine the distance to only one H2O maser source 
associated with SVS 13 in NGC 1333 in this paper, we 
should be able to reveal the overall 3-dimensional struc- 
ture of the Perseus molecular cloud complex through the 
further VLBI astrometry of other H2O maser sources in 
this region. 

5. Discussions 

5.1. Error sources in our astrometry 

The standard deviations of the post-fit residuals from 
the least-squares analysis, 0.10 mas and 0.29 mas in right 
ascension and declination, respectively, arc significantly 
larger than the formal errors in the Gaussian fitting of the 
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Fig. 4. Position measurements of the maser features associated with SVS 13. (a) The movement of the maser features in right 
ascension as a function of time, (b) The same as (a) in declination, (c) The movement of the maser features on the sky. Solid 
and dashed lines represent the best fit model with the annual parallax and linear proper motion for the maser features 1 and 2, 
respectively. Open circles and filled squares represent the observed positions of the maser features 1 and 2, respectively, and small 
crosses represent predicted positions of the maser features. Note that the formal errors given by the Gaussian fitting with the AIPS 
task JMFIT are as small as 0.1 mas or less, which cannot be shown in the figure. The reference position is the same as in Figure 3. 



maser features, 0.03-0.1 mas, implying that some system- 
atic errors significantly affect our astrometry. Although 
it is difficult to estimate the sources of the systematic er- 
rors in the VLBI astrometry quantitatively (Honma et al. 
2007; Hirota et al. 2007), they are mainly due to (1) the 
difference in the optical path lengths between the target 
and reference sources caused by the atmospheric zenith 
delay residual and/or (2) the variability of the structure 
of the maser feature for the following three reasons. 

First, the optical path length error due to an atmo- 
spheric zenith delay residual of 3 cm, which is typical for 
VERA observations (Honma et al. 2007), would cause a 
relative position error of 0.3 mas, for the case of SVS 13 
and J0336+3218 with a separation angle of 1.89° at an 
elevation angle of 20 degrees. This is comparable to the 
standard deviation of the post-fit residuals, especially in 
declination. Because the elevation angle of the sources 
was always higher than 20 degrees, the estimated position 
error of 0.3 mas gives an upper limit. Nevertheless, the 
atmospheric zenith delay residual could contribute signif- 
icantly to the error sources in our astrometry, especially 
in declination. 

Second, the internal structures in the maser features are 
found to be 0.04-0.27 mas based on our channel maps of 
the maser spots. Sometimes, the positions of the maser 
spots within a feature are spread over more than 1 mas, 
as found in the session 2005/021. Thus structural changes 
might explain the 1 mas declination residual in the data 
point of session 2004/355 in Figure 4. This data point is 
clearly an outlier. However, the magnitude of the possi- 
ble change in the feature positions due to the variation of 
the internal structure of the maser spots certainly could 
be comparable to the dispersion of the positions of the 
maser spots. This effect has also been found in the previ- 
ous results with VERA, in particular for the H2O maser 
sources associated with nearby molecular clouds (Hirota 
et al. 2007; Imai ct al. 2007). 

Third, the uncertainties in the station positions, delay 



model, and path length errors due to ionosphere would 
have negligible effects on our astrometry, according to the 
discussions in Honma et al. (2007). The uncertainties 
in the absolute position of J0336+3218 are reported to be 
0.46 mas and 0.56 mas in right ascension and declination, 
respectively (Fey et al. 2004). These uncertainties do 
not affect the derived annual parallax and proper motion 
because they would add only a constant offsets to the 
position of the maser features. 

In summary, we conclude that the astrometric accuracy 
in our observations is mainly limited by the atmospheric 
zenith delay error and/or the structure of the maser fea- 
tures. 

5.2. Origin of the H2O masers associated with SVS 13 

Along with the annual parallax, we successfully mea- 
sured the absolute positions and proper motions of the 
H2O maser features with VERA. As shown in Figure 5(b), 
the detected maser features are likely to be physically re- 
lated to VLA4A, only 0.2" (50 AU at the inferred distance 
of 235 pc) from the maser features, rather than VLA4B 
or other sources. The association of the maser features in 
the LSR velocity of 7-8 km s _1 with VLA4A is m agree- 
ment with the suggestion by Rodriguez et al. (2002). 
They also reported that another group of the H2O masers 
at the LSR velocities ranging from —15 to —25 km s _1 
is associated with VLA4B, suggesting the presence of an 
outflow from VLA4B (Rodriguez et al. 2002). Although 
the outflow activity of VLA4B cannot be excluded by our 
results without any detection of the H2O masers associ- 
ated with VLA4B, the present results show that VLA4A 
is the most favorable candidate for the powering source of 
the observed H2O masers. 

Figure 5(b) and Table 1 show the absolute proper mo- 
tions of the maser features associated with SVS 13. Both 
features are moving toward southeast with position an- 
gles of 114° and 133°, respectively, which are apparently 
in good agreement with the alignment of the HH ob- 
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Fig. 5. Distribution of H2O maser features associated with SVS 13 in NGC 1333. The reference position is the same as in 
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jects traced by the optical emission (Noriega-Crespo & 
Garnavich 2001) and molecular outflows traced by the 
CO J=2-l line (Bachiller ct al. 2000). 

Although the absolute proper motions of the maser fea- 
tures seem to be well aligned with the jets and outflows 
previously reported, we must stress that the proper mo- 
tions obtained with VERA do not represent the "intrin- 
sic" proper motions of the maser features because they are 
measured with respect to the Sun and, hence, include the 
contribution of the Solar motion. If we assume the Solar 
motion relative to the LSR based on the Hipparcos satel- 
lite data, ([/o,Vb,W )=(10.00, 5.25, 7.17) km s" 1 (Dehnen 
& Binney 1998), we can calculate the contribution of the 
Solar motion to the observed absolute proper motion to be 
3.7 mas yr -1 and —7.3 mas yr -1 in right ascension and 
declination, respectively. Subtracting these values from 
the observed proper motions listed in Table 1, the proper 
motions of the H2O maser features with respect to the 
LSR in right ascension and declination arc obtained to be 
(14.2, -0.6) mas yr -1 and (6.9, —2.7) mas yr" 1 , for the 
features 1 and 2, respectively, as shown in Figure 5(b). 
They are significantly different from those based on the 
absolute proper motions. 

Using the proper motions with respect to the LSR 
rather than the absolute ones, we will discuss the possible 
origin of the H2O masers associated with SVS 13. First, 
we rule out the possibility that Galactic rotation accounts 
for the observed proper motions, because the contribution 
of Galactic rotation would be only (—0.65,0.49) mas yr -1 , 
which is estimated by assuming R of 8.0 kpc (Reid 1993) 
and 9 of 236 km s" 1 (Reid & Brunthaler 2004). 

One of the most plausible explanations is that the 
proper motions of the maser features are due to the jet 



from VLA4A. In fact, VLA4A is proposed to be a power- 
ing source of the HH jets and molecular outflows as men- 
tioned above (e.g. Noriega-Crespo & Garnavich 2001; 
Bachiller et al. 2000). If VLA4A is at rest with re- 
spect to the LSR, the proper motions of the H2O masers 
with respect to the LSR, 15.8 km s _1 (14.2 mas yr -1 ) 
and 8.3 km s _1 (7.4 mas yr -1 ) for the features 1 and 
2, respectively, represent the motions relative to VLA4A. 
These observed transverse velocities are consistent with 
the jets and outflows previously observed in this region 
(Wootten et al. 2002). However, the directions of the 
proper motions arc not parallel to the jets from VLA4A 
(e.g. Noriega-Crespo & Garnavich 2001; Bachiller et al. 
2000; Wootten et al. 2002). In addition, positions of 
the maser features are not exactly aligned with the HH 
objects and molecular outflows from VLA4A. However, 
VLA4A itself might hava a proper motion with respect to 
the LSR, which could account for the apparent discrep- 
ancy found in the geometry of VLA4A and jets traced by 
the H2O masers. 

It is worth considering that the observed proper mo- 
tions are due to a rotating circumstellar disk of VLA4A. 
According to Anglada et al. (2004), VLA4B, a compan- 
ion of the close binary system including VLA4A, has a 
circumstellar disk with the size and mass of <30 AU and 
0.06Af Q , respectively, while the mass of the circumstellar 
disk of VLA4A is at least 5 times smaller than that of 
VLA4B. Note, however, that we detected a signature of 
acceleration of the maser features as shown in Figure 2, 
which might be indicative of rotating disk (e.g. Miyoshi et 
al. 1995). If we simply assume an edge-on disk with a ro- 
tation velocity v rot of 12 km s , which equals to the mean 
proper motion of the maser features, the radius of the disk, 
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r<i, can be derived from the relationship r c i = v^ ot /A to be 
18 AU, where A is the drift rate of the radial velocity 
of the maser features, 1.7 km s _1 yr" 1 . The inferred ra- 
dius of the disk is in good agreement with that of VLA4B 
(<30 AU; Anglada et al. 2004), but it is slightly smaller 
than the separation from VLA4A, ~50 AU. As a result, 
the enclosed mass within the radius of 18 AU is estimated 
to be 3M Q . This value gives lower limit of the mass be- 
cause we simply assume the inclination axis of the disk to 
be 0°. The mass of the envelope associated with SVS 13 
is derived to be of order IMq based on the intcrfcromctric 
millimeter wave observations of dust continuum emission 
(e.g. Looney et al. 2000). In addition, the luminosity of 
SVS 13 is as low as ~22i© (Jennings et al. 1987), sug- 
gesting that it should be a low-mass YSO. Although these 
results agree well with the mass of 3M Q estimated from 
the circumstellar disk model, the position of the disk (the 
maser features) and the YSO (radio continuum source) is 
inconsistent with the edge-on disk model. 

The results of our astrometry of the maser features are 
still puzzling because the positions and proper motions of 
the maser features with respect to VLA4A do not fully sat- 
isfy the requirements for either a jet or a rotating circum- 
stellar disk associated with VLA4A. One of the reasons 
for this discrepancy arises from the lack of the absolute 
proper motion measurements of VLA4A and other possi- 
ble candidates for the powering source of the masers. The 
proper motion vectors could be changed significantly by 
subtracting the absolute proper motion vector of VLA4A, 
making us reconsider all the possibilities discussed above. 
Other possibilities such as an orbital motion of the binary 
system consisting of VLA4A and VLA4B (e.g. Anglada 
et al. 2000, 2004; Rodriguez et al. 2002) could be 
tested if we measure the proper motions of VLA4A and 
VLA4B. Further high-precision astrometry of the maser 
features together with the possible candidate for their 
powering sources (i.e. VLA4A and other radio continuum 
sources in NGC1333) is essential to solve this problem 
(e.g. Rodriguez et al. 2005). 

6. Summary 

We present the results of multi-epoch VLBI astrometric 
observations of the 22 GHz H2O masers associated with 
SVS 13 in the NGC 1333 region with VERA. The principal 
results of this paper are summarized as follows. 

1 . We have determined the annual parallax of the H2 O 
masers associated with SVS 13 to be 4.25±0.32 mas, cor- 
responding to a distance of 235±18 pc from the Sun. 
Although the inferred distance of 235 pc is consistent with 
the photometric distance of 220 pc estimated by Cernis 
(1990), our results provide the distance to NGC 1333 with 
much higher accuracy than the photometric methods. 

2. The standard deviations of the post-fit residuals for 
the annual parallax and proper motion measurements are 
0.10 mas and 0.29 mas in right ascension and declina- 
tion, respectively. The astrometric error sources in our 
observations are discussed, and they are attributed to the 
difference in the optical path lengths between the target 



and reference sources caused by the atmospheric zenith 
delay errors and/or the variability of the structures of the 
maser features. 

3. The absolute positions and proper motions of the 
H2O maser features at the LSR velocities of 7-8 km s _1 
are derived, revealing that they are most likely associated 
with the radio continuum source VLA4A, based on the 
observation that the projected distances of the maser fea- 
tures are only 50 AU from VLA4A. 

4. We considered the possible origin of the observed 
proper motions of the maser features and found it difficult 
to explain the positions and the proper motions of the 
maser features in terms of either the jet or the rotating 
circumstellar disk associated with VLA4A solely based on 
the present results. Further highly precise astrometric 
observations of the maser features and radio continuum 
sources including VLA4A would be necessary to reveal the 
complex nature found in the HH 7-11 region in NGC 1333. 
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of Education, Culture, Sports, Science and Technology 
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